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We investigate the effect of quantum metric fluctuations on qubits that are gravitationally coupled 
to a background spacetime. In our first example, we study the propagation of a qubit in flat 
spacetime whose metric is subject to flat quantum fluctuations with a Gaussian spectrum. We find 
that these fluctuations cause two changes in the state of the qubit: they lead to a phase drift, as 
well as the expected exponential suppression (decoherence) of the off-diagonal terms in the density 
matrix. Secondly, we calculate the decoherence of a qubit in a circular orbit around a Schwarzschild 
black hole. The no-hair theorems suggest a quantum state for the metric in which the black hole's 
mass fluctuates with a thermal spectrum at the Hawking temperature. Again, we find that the 
orbiting qubit undergoes decoherence and a phase drift that both depend on the temperature of 
the black hole. Thirdly, we study the interaction of coherent and squeezed gravitational waves 
with a qubit in uniform motion. Finally, we investigate the decoherence of an accelerating qubit in 
Minkowski spacetime due to the Unruh effect. In this case decoherence is not due to fluctuations 
in the metric, but instead is caused by coupling (which we model with a standard Hamiltonian) 
between the qubit and the thermal cloud of Unruh particles bathing it. When the accelerating 
qubit is entangled with a stationary partner, the decoherence should induce a corresponding loss in 
teleportation fidelity. 

PACS numbers: 03.67.-a, 04.60.-m, 04.70.-s, 04.30.-w 
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I. INTRODUCTION 

For much of its history, experimental observation of 
quantum gravitational effects has been little more than 
an impossible dream for general relativity. Yet recent 
theoretical developments have begun to alter this view. 
New ideas in quantum gravity, including string theory, 
are being explored, which, if correct, would imply the 
presence of quantum gravitational behavior at length 
scales much larger than Planckian Lj L| Combined 
with emerging sensor and other technologies, these spec- 
ulations are raising the specter of near-term experimental 
tests for some of the most fundamental manifestations of 
quantum spacetime structure. 

An imperfect but illuminating analogy can be pursued 
in the physics of fluids. Consider a hypothetical stage in 
the development of the theory of fluids in which we are 
unaware of the atomic structure of matter, while we un- 
derstand the Navier-Stokes equations and are somehow 
aware of Planck's constant, h. In this analogy, the atomic 
structure of fluids is the unknown holy grail of quantum 
fluid mechanics which we are striving to discover with our 
limited knowledge. The fundamental quantities charac- 
terizing a fluid classically can be taken to be the speed 
of sound, c s , and the density, p. Along with H, these 
quantities can be combined to define a length scale 
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which would seem to set the quantum regime where clas- 
sical fluid mechanics must break down. For typical fluids 
(e.g. water), l q is of order 10 _8 cm; it does a remarkably 
good job of predicting the atomic scale. But in our hy- 
pothetical ignorance of atoms, we would be ill-advised to 
conclude that quantum effects must be negligible all the 
way down to the length scale l q . Indeed, we know (in 
hindsight) from kinetic theory that the true length scale 
which sets the breakdown of classical fluid mechanics 
is the correlation length, which is typically much larger 
than l q . Collective phenomena such as phase transitions 
may give rise to correlation lengths that are macroscopi- 
cally large. Even away from phase transitions, correlation 
lengths are large enough to have effects (such as Brow- 
nian motion) that are readily observable in experiments 
performed far above the lengthscale l q . 

What this analogy teaches us is to be open to the pos- 
sibility that, while the Planck length sets the scale for 
quantum gravity, there may be quantum effects of the 
unknown small-scale structure of spacetime which are 
analogous to phase transitions or Brownian motion, and 
which might be detectable at scales far above Planck. In 
this paper we will explore one possible source for such 
effects: quantum fluctuations in a background spacetime 
causing decoherence in qubits kinematically coupled to 
that background. In a sense, such decoherence would be 
analogous to Brownian motion caused by the small-scale 
structure of a fluid. We will make very few assumptions 
about the quantum theory of spacetime underlying the 
fluctuations, apart from demanding that it results in a 
Hilbert space structure for states of the metric which 
obeys the standard laws of quantum mechanics. We will 
conclude the paper by discussing a slightly different sce- 
nario for gravitational decoherence: the decoherence of 
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an accelerating qubit in Minkowski spacetime due to the 
Unruh effect. In this case decoherence is not due to fluc- 
tuations in the metric, but instead is caused by a cou- 
pling (which we model with a standard Hamiltonian) be- 
tween the qubit and the thermal cloud of Unruh particles 
bathing it. 

A separate set of motivations for studying gravitational 
decoherence arises from relativistic quantum information 
theory, a new field of research that studies the proper- 
ties of quantum information and quantum communica- 
tion as seen by observers in moving frames. One of the 
most important ingredients of quantum information the- 
ory is entanglement, and much of the interest so far has 
been directed to its properties under Lorentz transforma- 
tions. Czachor studied a version of the EPR experiment 
with relativistic particles and Peres et al. demon- 
strated that the spin of an electron is not covariant un- 
der Lorentz transformations Furthermore, the effect 
of Lorentz transformations on maximally spin-entangled 
Bell states in momentum eigenstates was studied by Als- 
ing and Milburn 6] , and Gingrich and Adami derived the 
general transformation rules for the spin-momentum en- 
tanglement of two qubits @ • Recently, these results were 
extended to the Lorentz transformation of polarization 
entanglement Q , and to situations where one observer is 
accelerated 0. Here, we take a look at a different aspect 
of relativistic quantum information theory: the effect of 
decoherence on a qubit due to quantum fluctuations in 
the metric. 

We will now turn our attention to the study of deco- 
herence in the four major paradigms of relativity: First, 
we will calculate the effect of flat fluctuations in the fab- 
ric of Minkowski spacetime on a linearly moving qubit. 
Next, we put a qubit in orbit around a black hole. The 
mass fluctuations due to the Hawking radiation induce 
fluctuations in the Schwarzschild metric, which in turn 
couple to our qubit. Following black holes, we will study 
the interaction of a linearly moving qubit with coher- 
ent and squeezed gravity waves. Coherent states corre- 
spond to classical gravitational waves which are expected 
to arise from astrophysical processes such as the inspiral 
of compact binaries, while squeezing arises due to the 
expansion of the universe. Finally, we study a related 
but qualitatively different decoherence phenomenon: we 
derive the decoherence of a linearly accelerating qubit in 
a bath of Unruh radiation. This treatment requires a 
field-theoretical description of the thermal Unruh bath, 
while we will be content with a first-quantized descrip- 
tion of the qubit, whose interaction with the quantum 
field of Unruh particles we will model with a standard 
Hamiltonian. 

First, we define a general quantum system to be in a 
state = J2k c k\k), where {|fc)}fc is a complete set of 
(non-degenerate and discrete) eigenstates of the Hamil- 
tonian. The free evolution of \ip) is given by 

\m) = E c k e iE ^ h \k) = £ c k e^\k) . (2) 

k k 



We consider the situation where the quantum system 
moves along a geodesic path in some fixed metric go, 
while the true metric g is subject to quantum fluctuations 
about go. The proper time of the system is denoted by 
r, and we have |'0( r )) = J2k Ck^ kT \k). The laboratory 
coordinates are t, x, y, and z. We want to know the state 
of the quantum system \ip(t)) in laboratory coordinates. 
This depends on the metric g: 

ds 2 — g^ v dx^dx" — dr 2 , (3) 

where the last equality follows from the geodesic motion 
of the quantum system. 

Now suppose that the metric itself is a (highly delo- 
calized) quantum system in a state p g . It is our central 
assumption that this quantum metric behaves as a reg- 
ular quantum system. In particular, we assume that the 
superposition principle is valid for quantum states of the 
metric. We initialize the state of our quantum system at 
time t = in \ip(0)) = J2k Ck \^)' ano - ^ ^ evolve freely 
within the metric pc for a certain period. The state of 
the system will generally become entangled with the met- 
ric, and we wish to determine the reduced density matrix 
of the system at time t. 

In general, the free evolution of an energy eigenstate 
in a particular quantum metric \gj) is given by 

\k)\ 9j ) ^ e^\k)\ 9j ) , (4) 

with dr = dt^J g^ v x v , and i M = dx^/dt. There is 
no back-action of the state of the qubit onto the met- 
ric. Specifying the (fluctuating) quantum metric and the 
geodesic motion of the quantum system allows us to de- 
termine the reduced density matrix of the quantum sys- 
tem. 

Rather than the general state in Eq. (J2J), we consider 
a qubit in a state 

|V) = -^(|0) + e^|l)) . (5) 

We could have chosen a more general state with differing 
relative amplitudes, but it turns out that this state gives 
rise to all the interesting physical features of the interac- 
tion model we study here. Similarly, the general case of 
an A^-level system does not give rise to conceptually new 
physics. 

Next, we will consider the geodesic motion of a qubit 
in three different fluctuating quantum spacetime met- 
rics. We consider "flat" fluctuations in a two-dimensional 
Minkowski space in Sec. [HJ mass fluctuations in the 
Schwarzschild metric for circular orbits of the qubit in 
Sec. IIIII and graviton-number fluctuations in coherent 
and squeezed gravity waves in Sec. IIVI 

II. FLUCTUATIONS IN FLAT SPACETIME 

In this section we calculate the density matrix of a 
qubit that experiences decoherence due to "flat" fluctu- 
ations in the Minkowski spacetime metric. 
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A. The Minkowski metric 

First, consider a uniformly moving qubit with velocity 
dx/dt = v in the ir-direction in a two-dimensional flat 
Minkowski space. Let the state of the metric be given by 



p g = daf(d) \gs)(9a\ 



(6) 



with a — (ai, a 2 , 03, 04) G K 4 and a particular metric 
given by 

_/ai+l a 2 W-l OUtti + l a 3 \ 

This defines the "flat" fluctuations in the two- 
dimensional Minkowski space. The distribution function 
f(a) determines the quantum fluctuations of the met- 
ric around a = 0. We assume that the fluctuations are 
Gaussian: 



/(a) = exp 
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ijAjkcik = exp [-(a, Aa)] , (8) 



where A is a positive symmetric variance matrix. The 
(dimensionless) parameters specified by the matrix ele- 
ments of A would ideally be derived from a complete 
quantum theory of gravity. In the absence of such a the- 
ory, we would treat the Ay as arbitrary real parameters. 

The proper time r of the qubit in motion, given a par- 
ticular metric gg, is given by 



where 



dr = dty/ g o + 2#ox v + g xx v 2 



9oo = a\ - aj + 2ai + 1 
gox = a,ia 3 - a 2 a A + a 3 - a 2 
gxx = 03 - a\ - 2a4 - 1 • 



(9) 



(10a) 
(10b) 
(10c) 



Furthermore, we assume linear motion, yielding dx — 
v dt, with H = c = G = 1. In the following, we will use 
Cl = u>i — u>q. 



B. The decoherence model 

The procedure to calculate the decoherence due to the 
above interaction of the qubit with the metric is as fol- 
lows: first we write the total state p = \tl)){ip\®p gi and we 
apply the interaction from Eq. (0J. We then substitute 
the expression from Eq. © into the resulting state. In 
fact we will use the polynomial expansion of the square 
root to second order. This is justified since the fluctu- 
ations are small. Subsequently, we trace out the metric 
state, since we have no direct access to its fluctuations. 
We can then evaluate the integral in Eq. 10, which yields 
an expression for the reduced density matrix of the qubit. 



The joint system of the qubit and the metric after the 
interaction of Eq. (0} is in the state 



|0)(0| + e- l(v+OT) |0)(l| 



-e 4(v+OT) |l)(0| + |l)(l| 



(11) 



Substituting Eq. into the above expression and ex- 
panding to second order will yield a state a|0)(0| + 
e~ llp P*\0)(l\ +e^/3|l)(0| +a|l)(l|, with 



= / dae- (s ^ = 



Vdet A ' 



(12) 



and 



(13) 

The proper time r(a) is approximated by 

r(o) =t[c + (c ll a) + (d,C 2 d) + O(d 3 )] , (14) 

with 

1 4 

— + . . . = Vl - v 2 = 7 (15a) 

(15b) 



co(v) 
ci(v) 



= 1 



V 

~2 



1-y ) (l,~v,v,-v 2 ) 

(i 



C 2 (v) 



2 

v 1 

2 
,2 



3^ v 2 
2 , 



V 2 \ 

, 3v 3 



(15c) 
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In order to calculate f3, we collect all the second-order 
terms a^a^ into the variance matrix A (yielding a new 
positive symmetric matrix B = A + i£ltC 2 ), and the 
first-order terms are collected in a linear exponential: 



(3 = e intl / dae- (S ' BS)+{il ' S) , 



(16) 



where u — ifM(l+v 2 /2)(l, —v, v, —v 2 ). The overall phase 
factor originates from the time dilation observed for a 
moving body with velocity v. The integral in Eq. 1|16|) 
can be evaluated formally to yield 



dSe -(S,BS) + (uM) = 



VdetB 



(17) 



For the integral to exist, the real part of the eigenvalues 
of B must be strictly positive, which is ensured by A. 

Given a particular variance matrix A, we can calculate 
the (normalized) state of the qubit as a function of the 
travelled coordinate time t: 



1 



1 



vj e 



1 



(18) 
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Re 



FIG. 1: The amplitude and phase of the off-diagonal ele- 
ment of the density matrix. The inward spiralling is due to 
the phase drift and the decoherence. The knee in the figure 
around the point (Re, Im) = (0.8,0.1) is due to the complex 
factor r\. 



where ip' = ip + 7f2f is the time-dilated free evolution 
of the qubit, and r\ = y/det A/ det B. Furthermore, we 
defined 



V 2 = Re[-(u,B~ i u)] and 5 = Im[- (u, B' 



Whcn the variance matrix is diagonal (^4 = a~ 1), i.e., 
when there are no special correlations between the differ- 
ent space-time dimensions, we can give a fairly straight- 
forward expression of the decohered qubit. Since the 
quantum fluctuations of the metric are assumed small 
(Planck scale), the variance a is small, and we can ex- 
pand the solution around a = 0. Using Mathematica, we 
find that 



r = (i + v 2 )nta 



(20) 



and 



S = -v 2 (5 + 5v 2 + llv 4 + 3u 6 ) n 3 t 3 a 4 . (21) 

Interestingly, there are three effects that alter the state 
of the qubit. First, there is the expected decoherence 
exp(— r 2 ), the argument of which scales quadratically 
with time t, frequency difference O, and variance a. Sec- 
ondly, we found a phase drift 5, due to the interaction 
with the quantum metric. This effect scales with (fit) 3 
and a , and it arises when there are off-diagonal terms 
in the matrix B, or, equivalently, when there are ajdj- 
terms with j ^ k in Eq. 11411 . Finally the factor rj be- 
haves according to 



V 



^4 « 1 + ^(1 + v 2 ) 2 Vlta 2 
det B 2 y ' 



(22) 



The evolution of the complex off-diagonal element is 
shown in Fig. ^ 

The interaction of the quantum fluctuations of the met- 
ric with a qubit generalizes directly to the case of an 
iV-level system, or quNit. The density matrix will be an 
NxN matrix, and the off-diagonal elements will be of the 



formrye^fc+^e r j fc . We obtain Lp'- k , Sjk, and F 
Eqs. {TBJ|, 129 and (|2U by substituting Qjk 

for n. 



from 



C. Ensemble of qubits 

Note that the above calculation describes an experi- 
ment in which the ensemble of qubits necessary to statis- 
tically measure the reduced quantum state of the qubit 
is obtained via repetition of each ensemble qubit 's time 
evolution through the same metric state Eq. |SJ What 
happens when we send an ensemble of N qubits with 
uniform velocity v, while their states are all subjected to 
the same fluctuations in the Minkowski metric simulta- 
neously? Let \s) be an JV-bit string, and < n(s) < N 
the number of ones in \s). The evolution of a bit string 
\s) on the metric \gs) in terms of the proper time r is 
given by 

\s)\gs) -> exp[z?i(s) flT S }\s)\g s ) ■ (23) 
There are 2 N such strings, comprising a set iS, and the 



(19) density operator becomes a 2 x 2 matrix: 



Pn 



dae 



-(a, AS) 



E 



J[n(s)-n(s')nt] 



s)(s'\ 



(24) 



It is immediately obvious that the states most sensitive 
to this type of decoherence is the maximally entangled 
GHZ state |00 . . . 0) + 1 1 1 . . . 1), that is, states that have 
maximal n(s) — n(s'). As a direct consequence, these 
states can in principle be used to detect the fluctuations 

Similarly, it is now also clear that such decoherence 
can be countered by encoding quantum information in 
the decoherence-free subspace spanned by the subsets of 
S for which n(s) — n(s') — 0. For a given n(s), there are 
( , such states. This behaviour is reminiscent of the 
technique used in Ref. 11]. 



III. MASS FLUCTUATIONS IN BLACK HOLES 
A. The Schwarzschild metric 

Now, we will turn our attention to the situation of a 
qubit orbiting a black hole, and we again assume that 
the metric is a quantum mechanical object. The no-hair 
theorems for black-hole structure suggest that quantum 
fluctuations in the metric will consist of fluctuations in 
the mass of the black hole (hence we neglect the fluc- 
tuations in angular momentum, and assume that charge 
fluctuations are forbidden in quantum gravity because of 
super-selection rules) . 

For simplicity, we assume that the orbit of the qubit is 
circular. The Schwarzschild metric of a given mass M is 
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given by 



ds 2 = g M = - 1 - 



2M 



dt 2 



dr 2 



1 - 



2M 



2 dn 2 , (25) 



where dVL 2 = d9 2 + sin 2 9 d<j) 2 . 

The evolution of energy eigenstates is again given by 
Eq. 10}. We need to express r in terms of t. To this end, 
we start from the Hamiltonian for geodesic motion in the 
Schwarzschild background 



H 



2 _ AM + 

r 

E 2 

- 2 (1_2M) 
1 



2M 
r 

1 - 



p 2 . 

2 
2M 

r 



2r 2 



2r 2 sin 2 9 
2r 2 



(26) 



Here we have assumed without loss of generality that the 
orbital plane has 9 = ir/2 and thus by spherical sym- 
metry pq = 0. Furthermore, we introduced the energy 
E = —po and the angular momentum L = p$. The last 
equality holds since the four- velocity v = dx/dr is nor- 



malized: 



-1. Furthermore, we have 



Pa = 9ooV 



2M\ dt 
d7 



-E , 



(27) 



which gives us an expression for r in terms of t and E. 
We now need to determine E. 

^From Eq. (|26l) and p r = (since we consider a circular 
orbit), we derive 



E 2 



r j 



V(r) 



(28) 



For stable orbits we need the potential V(r) to be a min- 



dV_ 

dr 



and 



3 2 V 



dr 2 



> 



This yields 



( 2M\ L 2 M / L 2 \ _ () 



Mr 2 
r - 3Af 



(29) 



(30) 



(31) 



which makes manifest the well-known fact that there are 
no stable circular orbits inside the critical radius r = 3M. 
For any ro > 3M, a circular orbit at r = ro can be 
found with the above value of L. On such an orbit, the 



energy E is determined by substituting L 2 into V(ro) 
from Eq. (27) we find that 



dr 2 = 1 - 



3M 



dr 



(32) 



We will assume that the black hole is in thermal equi- 
librium with the thermal Hawking radiation it emits at 
the temperature Th- This is equivalent to assuming that 
the black hole is inside a conducting sphere with per- 
fectly reflecting walls placed at a large radius away from 
the horizon. The quantum state of the metric can then 
be written as in Eq. JHJ in the form 



dM exp 



(M — M ) 2 
2a 2 , 



\9m){9. 



M 



(33) 



where M denotes the classical mass of the hole, about 
which the mass M fluctuates in equilibrium with the 
thermal bath at temperature Th- The fluctuations are 
Gaussian (in accordance with the thermodynamic limit) 
around Mq. More precisely, a general canonical ensemble 
with Hamiltonian H has energy fluctuations given by 



(H 2 



(H) 



dp- 



log Z 



d_ 

dp 



(H) , (34) 



where Z = J e~^ H is the partition function, j3 — l/(fcsT) 
is the inverse temperature, and (• • • ) denotes the canoni- 
cal ensemble average; thus, e.g., (H) is equal to the inter- 
nal energy U. According to our assumption as explained 
below, afj will be given by 



2 



(H 2 



'Th 



(H) 



(35) 



i.e. equal to the energy fluctuations of a Bose gas at ther- 
mal equilibrium with black-body radiation at the Hawk- 
ing temperature. 

In order to find the effect of thermal fluctuations in 
the black-hole mass on the evolution of coherence in the 
state of our qubit in circular orbit at r — ro, we calculate 
the off-diagonal element of the qubit 's density matrix: 



dM exp 



9M 2 flt 



8r 2 + m{VL) 2 t(j 2 M 



(M-Mo) 2 ] tflt,/l-2M 

—^5 e * r ° 

exv(-T 2 + iS) , (36) 



where 



_ 9s/2 /3M + 2r \ ... 

1 = — — I ^ I SiCCTM 



(37) 



and 



S = QU 1- 



3M 



81 {3M + 2r x o3+3 a 



128 V rg 



VL 6 t A a%, . (38) 
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The first term of the phase drift is the regular relativis- 
tic effect, and the second term is the variational phase 
drift due to the quantum fluctuations in the mass of the 
black hole. The off-diagonal element of the density ma- 
trix evolves similar to the one in the previous section (see 
Fig.QJ. 



B. Mass fluctuations 

Now we seek an expression for um . We know that the 
Hawking temperature of a black hole is proportional to 
the surface gravity k [l3j : 



kT „ = h K = Hc 

c2n 8nGM ' 



(39) 



where we have used the fact that a non-rotating neutral 
black hole of mass M has K = (AGM/c 4 )^ 1 , and therefore 
I3- 1 = kT H = hc 3 /(8nGM) \rft. For clarity we will 
keep the physical constants H, c and G in arbitrary units 
throughout this subsection. 

The Hawking radiation has a perfect black-body spec- 
trum. We will assume that the fluctuations in the mass of 
the hole are identical to the fluctuations in the average 
energy U of a black body of surface area 4ttG 2 Mo 2 /c 4 
(the area of the black-hole horizon) radiating at the 
Hawking temperature Th- The average internal energy 
of a black body (Boson gas) of volume V at temperature 
T is given by 



(H) = U = 



4er 7r 2 /c B 4 
— VT = — 

c 15h 3 c 3 



VT 4 = 



V 



15fi 3 c 3 /3 4 



(40) 



where a is the Stefan-Boltzmann constant. According to 
Eqs. (34)-(35) 



4 2 



OP 



{Hf 



4tt 2 V 
I5h 3 c 3 ~W 



(41) 



Note that the black-body energy U itself is not equal to 
the classical hole mass Mo; it is only the mean-square 
fluctuations a 2 M about M that we assume are identical 
to the thermal fluctuations in U at the Hawking tempera- 
ture. Substituting the Hawking value (3 = 8irGM^j '(he 3 ), 
and the geometric identity V = (32vr/3)G 3 M 7c 6 in 
Eq. (41), we finally obtain 



To give a feeling for the order of magnitude of this 
decoherence effect, we note that in order to have the 
off-diagonal elements of the qubit state decohere to (an 
absolute value of) 1/e, the value of T [cf. Eq. (37)] 
must be 1. Given a black hole of one solar mass (10 38 
Planck masses), and a qubit with transition frequency 
O = 10 15 Hz in the lowest stable orbit (r = 3M ), the 
(asymptotic, Minkowskian) coordinate time necessary to 
reach this level of decoherence is 10 13 years. To get a 
decoherence effect of the same magnitude over a time 
interval of one year, the black hole must be no more 
massive than 10 kilograms. On the other hand, if we 
take the suggestions of microscopic black-hole produc- 
tion in the hadron colliders currently under construction 
seriously [l| , their decoherence effect on nearby quantum 
states should be considerably larger. 



IV. GRAVITATIONAL WAVES 

In this section, we will consider the interaction of a 
qubit with coherent and squeezed graviton fields through 
the proper time evolution of the qubit. In the "transverse- 
traceless" gauge, the metric for classical gravitational 
waves travelling in the z-direction is given by |15| 

ds 2 = -dudw + (l + h c + ) 2 dx 2 + {l-h c + ) 2 dy 2 

+2h c x dxdy, (44) 

where u = t — z and w = t + z, and h c + x the classical 
amplitudes of the gravitational wave in the + and the 
x polarization. Furthermore, the qubit moves in the xy 
plane with velocity v such that 



du dw 
~d±~ ~dl 



1 



dx 
~dl 



v sm( 



dy 
dt 



vcos9 . (45) 



This defines 9. Expanding the metric to linear order in 
h° + , the proper time of the qubit is then given by (dr 2 = 
ds 2 ): 



T = tJ(l- v 2 ) - 2h c + v 2 cos 26* -I- h L x v 2 sin26>. (46) 

When we quantize the (weak) gravitational wave, the 
metric can be approximated by 



2 

a M = 



1 



11520 tt 2 G 2 M n 2 11520 it 2 V M< 



(42) 



ds 2 = -dudw + (1 + 2nh + )dx 2 + (1 - 2nh + )dy 2 

+2mh x dxdy, (47) 



where m p = y/hc/G is the Planck mass. Going back to 
geometric (Planck) units (where Planck mass is unity), 
we can write Eq. (42) in the form 



where n,m are the graviton numbers for h + and h x re- 
spectively. The classical metric is then retrieved by set- 
ting 



1 



' M 



1 



11520 7T 2 VM| 



(43) 



which is the value that needs to be substituted in 
Eqs. (33)-(38). 



h° + = (n)h + and hZ = (m)h x 



(48) 



Here (n) and (m) are the average graviton numbers. The 
h + and h x denote a distribution functions over a set of 
frequencies, and may be a function of u and w. 
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The proper time of the qubit interacting with n +- 
polarized and m x -polarized gravitons, expressed in the 
standard coordinates, then becomes 

dr nm = dt\J (1 — v 2 ) — 2nh + v 2 cos 29 + mh x v 2 sin 29 . 

(49) 

Suppose that we deal only with plane waves of specific 
frequencies lo + and u> x such that h + (t,z) — /i+cosw+t 
and h x {t 1 z) = h x cos(uj x t + ip). We can then approxi- 
mate this expression using \J a + b « \fa + b/(2y/a): 

ft 

T~nm = — n / dt' 2 h+v 2 cos 29 cos uj + t' 



+m dt' h x v 2 sin 29 cos(u> x t' + <p) . (50) 
Jo 



with 7 = \/\ ~ v 2 . We can compactify this expression 
and write r„ m = jt — nr + + mr x , where 



2h + v 2 cos 29 sin oj + t 
h x v 2 sin29 sin(tu x t + tp) 

JUJ X 



Let's consider quantized gravitational waves in the two 
polarization modes h + and h x . That is, we can write 
the state of the mode in terms of bosonic creation and 
annihilation operators h, and hj, where j G {+, x}, and 



[hj,hl]=5 jk and [h j ,h k ] = [h],hl]=Q. (52) 

We will now use this formalism to study the effect of 
different quantum states of gravitational radiation on the 
qubit. 



A. Coherent gravitational waves 

The most classical states of a bosonic field are the co- 
herent states, so it is natural to study coherent gravi- 
tational waves 0| . We are looking at single modes of a 
single frequency. Analogous to optical coherent states we 
define the coherent state of the (polarized) gravitational 
wave with amplitudes rjj (J G {+, x}) as 



Wi) 



00 v? 



where |n) is the n-graviton Fock state obtained by 



hj\n}j = \fn + I|n + l)j and hj\n)j = y/ri\n — i)j . 

(54) 

Here, we have chosen rj + , rj x real without loss of general- 
ity. 

When a qubit interacts with the gravitational wave, 
we write the state of the joint system as 

|0) + |1) ^ | )+e'"^|l) 



V2 



n,m— 



V2 



V n\m\ 



(55) 



Tracing out the gravitational wave then leaves us with a 
qubit in the state 



1 (3 
f3* 1 



where 



/3oc 



°° „2n -inOT| 



n.m— 



(56) 



(57) 



r] x }. This yields 



up to a constant exp[ijQt — i 
/3 = e^ nt exp [T 1 2 + {e- %nT + -1) + T)l (e^* ~ 1)] ■ (58) 

Expanding the the phase factor according to e %x = cos x+ 
isinx, we can formally rewrite (3 = e~ r+lS , with 



(51a) 



(51b) and 



T = rf + (1 - cos</>+) (1 - cos^x) 



vl4>l 



5 = Qt'f — rn_ sin <j) + + rj 2 , sm < 

2 l ,2 



Here, we defined 



and 



2h+v 2 cos 29 sinu! + t 
uj + 7 
Q, h x v 2 sin29 siri(uj x t + ip) 

lu x 7 



(59) 

(60) 

(61a) 
(61b) 



This means that there is a periodic decoherence of the 
qubit under influence of the gravitational wave. Indeed, 
the restoration of coherence is induced mathematically 
by the occurrence of a phase factor in the argument of 
the exponential. Similarly, the induced phase drift S is 
periodic in time. However, since <f>+ and <p x are them- 
selves periodic functions of the laboratory time t with 
a small amplitude [see Eq. (|51|) ]. the periodicity of the 
decoherence and the phase drift will be extremely diffi- 
cult to observe in practice, due to their extremely small 
amplitude. 

However, in practice the qubit will not be interacting 
with a plane wave, but with a pulse of finite extension. 
If we consider a simple coherent Gaussian wave packet, 
the gravitational wave in one polarization can be written 
as 



|*(t)> = —_£ 



duj i 



(62) 

For simplicity we have included the dependence of r/ on 
lu in the Gaussian function. 

Under the evolution \k) <X> \n) g — ► e inkTn \k) <x> \n) g the 
metric and the qubit state \ip) = (|0) + |l))/\/2 become 
entangled. When we trace out the gravitational wave, 
the state of the qubit then becomes 
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Pout 



V47TCT 2 



/ <iw da;' exp 



2ct 2 



2ct 2 



— u> )t 



Vina 2 



n,m— 

+ e- inT "^|0)(l| +e mr ^\l)(0\ + |1)(1| 



E/ Q ^^ 1 ^^(|0)<0|+ e -^"^|0)<l|+ e ^^|l)<0| + |l)<l|) = ^ , (63) 



\Jn\ ml 



\n u ){m u 



a (3* 



where r n (cj) = nTk(uj) with given by Eq. (|51|l and 
lu = ujk- We can calculate the matrix elements of p ou t : 




with jit = 2hv 2 cos 29 /^v. In order to calculate /?, we 
note that n/i is of the order of the classical amplitude 
of the gravitational wave, i.e., 10" 17 at best. We can 
therefore approximate the imaginary part of the expo- 
nential to the first few orders, using e x ~ 1 + a; + x 2 /2 
(x <C 1). We need the second order in the approximation 
of the exponential, because it is going to contribute to 
the decoherence [Tr(p 2 ut )] to the same order as the first 
term: 



a 1 . irf^ 
P - 77 



2 \J Ana* Jo 



du> — smut e ^ 



77V 2 r °° 



2V47Tt7' Jo 



auj I — sin tot 



(<->—> r 



(66) 



The expression rj 2 /j, is on the order of the classical gravitational wave amplitude h c , with rj ^> 1 and fi <gC 1. Now, /3 
can be written as 



1 I ih c n^ ,/(rf\ /i c2 fi 2 V^ 



2rr 



erf 



2(7 



'7rcr 



+ f erf (at) 



(67) 



The imaginary part of j3 is shown in Fig. [21 The magnitude of the decoherence [given by Tr(p 2 ut )] can now be 
calculated: 



Tr(p 2 ut ) = 1 - 



h c2 n 2 ir 
~2a 2 ~ 



: erf(cri) erf 

'it 2 



(68) 



The general behaviour of the decoherence is shown in Fig. [21 



B. Squeezed gravitational waves 

Now let's look at squeezed gravity waves. These states 
are expected to arise due the expansion of the uni- 
verse |l8j . Suppose the squeezed gravitational wave in 
one polarization mode (for simplicity) has the following 
graviton-number expansion: 

71=0 ISI 

where £ is the complex squeezing parameter, which de- 
pends on the rate of expansion of the universe. 



A squeezed gravitational wave packet can be defined 
as (tanh |£| = r): 



IC(*)) = 




(70) 

Using the techniques of the previous section, we can ex- 
press the state of the qubit again in the form of 

HG 1 *?)' (71) 
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where 
= 



1-r 2 f°° e-^^^duj 

Ana 2 J y/\ — e i4>(^) r 2 



r 2 (2 + r 2 ) 
X 8(l-r 2 ) 2 a 2 



rf 



at 



1 + V7rcrt erf 



at 



(72) 



where we defined 4>{oj) — /iQ sin ut/uj. The off-diagonal 
elements of the density matrix therefore exhibit the same 
behaviour as in the coherent case. The field amplitude of 
the squeezed gravitational waves are now proportional to 
[ir 1 j{\ — r 2 ), where < r < 1 is the squeezing parameter. 
The decoherence with time is again given in Fig. [21 



LFtW 




RRW 



FIG. 3: The world line of an accelerated observer. LRW and 
RRW denote the left- and right-Rindler wedges, respectively. 
They are the spacelike separated regions with respect to the 
origin. The accelerated observer intersects the x-sods at xq. 



V. DECOHERENCE VIA UNRUH RADIATION acceleration along the worldline (see Fig.©: 



This section is different from the others in that we will 
not consider a fluctuating quantum metric, but rather 
an accelerated qubit in Minkowski spacetime interact- 
ing with a thermal bath of Unruh particles. The quan- 
tum field whose states contain the Unruh particles and 
which interacts with the qubit will be assumed to be a 
real scalar field for simplicity, but qualitatively similar 
results are valid for any linear field independently of spin 
or charge. When the qubit is an atomic system, the rel- 
evant quantum field is that of photons in flat spacetime. 



A. Acceleration in flat space-time 

We closely follow the treatment given in ^3 °f t ne be- 
havior of an accelerating particle detector in Minkowski 
spacetime. Accordingly, we consider a qubit in uniform 



g 0.95 

c 

<D 
M 
<D 

o 
o 

0) 

a 



. 9 



. 85 



. 8 




Time 



FIG. 2: The behaviour of Tr(pSut) for different a. When a 
becomes larger, the trace remains closer to one. The values 
of a are 1, 2, 3, and 5. 



t = — sinh(ar) and x = — cosh(ar) , (73) 
a a 

where r denotes proper time along the worldline, and a 
is the magnitude of the acceleration. As in all standard 
treatments of the Unruh effect, we envision a congruence 
of uniformly accelerating world lines, where the world- 
line crossing the a;- axis at x = xq has acceleration 1/xq. 
The coordinate origin can be adjusted so that our given 
qubit 's worldline crosses the x-axis at x = 1/a in ac- 
cordance with Eq. J73J). We will assume that the qubit 
interacts with a real, massless scalar field <p propagating 
on Minkowski spacetime in just the same way as a "par- 
ticle detector" for <fr would. Accordingly, we will assume 
a model interaction Hamiltonian 



Hi(t)=e{t) / 4>{x,t) 4>{x)b + 4>*{x)P 



-g (fx . 

(74) 

Here the integration is over the global spacelike Cauchy 
surface E = {t = constant} in Minkowski spacetime, 
and e(t) is a coupling constant which is explicitly time 
dependent to allow the acceleration to vanish outside a 
finite time interval; we will assume that e(t) is a constant, 
e(t) = e, within a finite interval in t of duration A and 
zero outside that interval. The function ip(x) is part of 
the coupling constant, ft is a smooth function that van- 
ishes outside a small volume around the qubit, and that 
models the finite spatial range of the interaction. 

The field operator 4>(x, t) can be expanded either as 
the mode sum 

$(x,t) = ^2 \ui{x,t)a M {ui) +Ui(x,t)*a\, 1 (u 2 ) , (75) 

i 

where {ui} is an orthonormal basis of positive-frequency 
(with respect to the timelike Killing field d/dt) solutions 
of the (Klein-Gordon) field equation for <j> on Minkowski 

spacetime and a^ M {ui) and a,M(ui) denote the corre- 
sponding creation and annihilation operators, or as the 
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mode sum (H.c. is the Hermitian conjugate) 



B. Evolution of the qubit 



t) 



[w m (x,t) a R (w Ri ) + H.c.] 



(76) 



in the right Rindler wedge {x > \t\} of Minkowski space- 
time (where the accelerating qubit 's worldline is con- 
tained), where {w R i} are an orthonormal basis of so- 
lutions on the right Rindler wedge which are positive- 
frequency with respect to the timelike boost Killing field 
a{xd t + td x ) there (which is normalized to have unit 
length along our qubit 's worldline) and a R (w R i) and 
a R (w R i) denote the corresponding creation and annihi- 
lation operators. 

There is a similar mode sum for <j)(x, t) on the left 
Rindler wedge (see 0, f°r details) . The operators V* 
and b in Eq. 1741) denote the raising and lowering oper- 
ators in the internal (two-dimensional) Hilbert space of 
the qubit: 



6|0) = 0, 6|1) = |0) 
St|0) = |l), &t|l>=0. 



(77) 



Accordingly, the internal Hamiltonian of the qubit is sim- 
ply 



H = nPb , 



(78) 



where is the energy difference between the ground and 
excited states. 

As first discovered by Unruh [14|. when the (bosonic) 
Fock space for the quantum theory of the field <f> is built 
from the Rindler mode expansion Eq. I)76|) instead of the 
Minkowski- mode expansion Eq. I|75l) , the Minkowski vac- 
uum state |0)m can be written in the form 



10) 



M 



(79) 



Here the product is over all Rindler modes of the form 
uum (and the corresponding left-wedge modes wli) with 
(positive) frequencies tOi, and the inner sum is over all 
non-negative integers rii with \rii,R} and \rii,L) denot- 
ing the state with rii particles in the mode wm and wn 
respectively. The normalization constants Cj are given 
by 



(80) 



When the left-Rindler components of the vacuum state 
|0)j\/ of Eq. H79fl are traced-out, the reduced density ma- 
trix for Minkowski vacuum as seen by uniformly accel- 
erating observers in the right Rindler wedge is given by 

PR = II (°* 2 E e~ 2mWa K R) ® (ni, R^j (81) 

which is exactly a thermal state at the Unruh tempera- 
ture 



k B T = — 
2ir 



(82) 



We are interested in calculating the final state of the 
qubit when initially, before the acceleration commences, 
it is in the state 



10) 



ID 



V2 



(83) 



Let us first assume that initially the qubit is in the 
ground state |0), and the field is in the Fock state |n) 
having n particles in the mode corresponding to a one- 
particle Hilbert space element (positive-frequency solu- 
tion) x(x, t). So the combined initial state of the field 
and the qubit is 



|in) = |n) ® |0) 



(84) 



The evolution of this input state under the interaction 
Hamiltonian Eq. I|74|) is computed in first order pertur- 
bation theory (in the coupling constant e) in [l2] | , between 
Eqs. (3.4) and (3.26) there. To state their result, let us 
introduce the smooth function of compact support 



(85) 



and the retarded minus advanced solution of the Klein- 
Gordon equation with source /: £f = Rf — Af, where 
□ (i?/) = □(A/0 = /• Denote the positive and negative 
frequency parts of Ef by T + and T_, respectively, as 
in [13] ■ We have the basic relation 

kf) = J mt)f(x!,t)^d 4 x 

= a(T*_ ) - fit(r+) (86) 

expressing the field operator smeared with the test func- 
tion / in terms of the creation and annihilation operators 
acting on Fock space. Then, at the end of the acceler- 
ation (when t > A), the final state of the field and the 
qubit system evolving under the interaction Hamiltonian 
Eq. Ij74(l is, according to first-order perturbation theory, 



|out) = \n) <g) |0) -ia(T*_)\n) ® |1) 



(87) 



Here 



a(r)|n) - v^(r*,x)|n-i) 

= Vn f(x 1 t)x(x,t)^/^gd 4: x\n- 1) 



where (, ) denotes the inner product on the 1-particle 
Hilbert space of positive-frequency solutions. Assume 
now that the initial state is 



\n) (g) |1) 



(89) 



instead of Eq. I|84|l . The evolution of this input state can 
be computed in exactly the same way as in 12], with 
only slight modifications of their Eqs. (3.18) to (3.25). It 
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follows that at the end of the acceleration (when t ^> A), 
the final state evolving from Eq. I|89ll according to first- 
order perturbation theory is 



|out) = \n) ® |1) + ia){T+)\n) ® |0) , (90) 



where 



a t (r+)|n) - yfc+lv |f + <E> S (n • x)) , 

v = v/(r+,r + ) , f + = r+/i/. (91) 

Here |I + (E>s (n • x)) denotes the Fock state that cor- 
responds to the symmetrized product of the n + 1 nor- 
malized positive- frequency solutions (one-particle states): 
r + • x ■ X ' ' ' X- Putting together Eqs. I|87|l and (|90|l . when 
the initial state of the field and the qubit system is the 
coherent superposition 



|in) = |n)®(|0> + |l»/V2, 



(92) 



the outgoing state at the end of the accelera- 
tion/interaction episode is 



| out) 



|n)®(|0) + |l»/V2 

-*Vn(T* )X )|n-l)®|l|)/V2 

+iv Vn + T | f + ®s (n • x) ) ® |0)/V2 .(93) 



The physical interpretation of Eq. (|93(l is simple: the first 
term is the zeroth-order unperturbed initial state, the 
second term corresponds to the absorption of a single 
quantum of ^-radiation by the qubit, and the third term 
corresponds to the stimulated emission of a single quan- 
tum at the transition frequency f2. When the initial state 
of the field is a more general Fock state containing n\ 
particles in mode xi , n 2 particles in mode xi an d so on: 



Eq. I|76(l . which is the natural construction of the quan- 
tum field theory of <p from the viewpoint of an observer 
accelerating with the qubit along the same worldline. Ac- 
cording to such an observer, the Minkowski vacuum state 
in the right Rindler wedge is given by the thermal density 
matrix Eq. (|81|) . The initial state of the combined field 
and qubit system at the start of the acceleration is 



Pin = PR® \ip){ip\ 

= uPR® (|0) + |1))®' 



+ <1|). (96) 



We will now make the formal simplification that the 
thermal state Eq. I|81(l contains only one Rindler mode 
w R.i = Xfi? where xn is a mode whose frequency to equals 
the transition frequency Q. This simplification is justi- 
fied because the overlap factor (r*_, x) in the absorption 
(second) term in Eq. (|93|l is negligible unless x happens 
to be a mode at a frequency to w f2 [see Eq. 1(88(1 ] . While 
the spontaneous emission terms into states with other 
modes are not similarly negligible, they have a trivial 
form whose contribution is qualitatively the same with 
or without the simplification. With this assumption, the 
product over modes in Eq. I|81|) disappears, and the ther- 
mal state of Rindler particles bathing our qubit can be 
written in the form 

PR = C 2 Y, e- 2 ™°/>, R) (n, R\ , (97) 

n 

where 



c= vT 



-27rO/a 



(98) 



The input state of the field and the qubit system, 
Eq. I|96(l . now takes the simpler form 



|m) = K,n 2 ,--- ,n,)®(|0) + |l))/^, (94) 

it is not difficult to calculate the appropriate generaliza- 
tion of Eq. (|9"5|l for the outgoing state; the result is: 



(out) = |ni, n 2l 
q 



^^\ / ^( r -'X fe )l ni '"' 

v 1 fe=l 



n q ) (|0) + |1))/V2 

,n,)®|l) 



V2\ 



fc=i 



n k |f + ®s(m • xi) 



(» 9 -X 9 ))®|0) ■ 
(95) 



We now have all the machinery we need to compute the 
evolution of the internal state of an accelerating qubit 
which is initially in the state Eq. I|83|l . Notice that at 
no point in the above formalism, Eqs. I|84() to 1)95(1 ■ we 
made any reference to the particular construction of the 
Fock space for the field <fr. In particular, the formalism is 
valid for the Rindler construction based on the mode sum 



C 



= T ^e- 2 ^[K J R)®(|0) + |l))] 

n 

®[<n,iJ|®((0| + <l|)] . 



(99) 



Each of the components in the sum in Eq. I|98|l undergoes 
an evolution as in Eq. I|95|) . ending up at the end of the 
acceleration in the final state: 



|0) + |1) 



1 



In, R) ^ 

tVn(T*,xn)|n-l,iZ>|l) 



|n,i2)(|0) + |l)) 



+ %vy/n + l\T + ® s (n • xo) , R) ® |0) , (100) 



where we used 



l + -\{V^Xn)\ 2 + {n+l)- 



(101) 



We have corrected the normalization error affecting 
Eqs. (|87|l through l|95|l which is symptomatic of any ap- 
proximation based on first-order perturbation theory. 
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FIG. 4: The purity Tr(p| out ) as a function of the acceleration 
of the qubit. The acceleration is plotted in units of Q (i.e. the 
horizontal axis depicts the dimensionless quantity a/f2).When 
there is no acceleration the qubit remains pure (purity= 1), 
and for increasing acceleration the qubit approaches the max- 
imally mixed state (purity= 0.5). 



The crucial quantity measuring the strength of the in- 
teraction is the overlap fi = (T*_,xn)> which, according 
to Eqs. 1(85(1 and 1(88(1 . can be calculated as 



V2n 

eA 



<P*(x)e ik °- s d 3 x 



20 



(102) 



where ko denotes the wave vector corresponding to the 
Rindler mode \(i (thus |fco| ~ ^, at least for small accel- 
erations a), k is a length scale setting the spatial range 
of the interaction (thus ip(x) falls off like a Gaussian dis- 
tribution with variance k 2 ), and ^ is a phase encoding 
directional information about ko which is irrelevant to 
our considerations. Similarly to Eq. 1(102(1 . the quantity 
v [Eq. (frlTJl ] can be computed as 

(103) 



2(v^3) 



1/2 



Getting back to the evolution of the input state 
Eq. |J5SJ|, combining Eq. (|100|l with Eq. gives 



Pout 



C 2 Er 



— 2-TrnSl/ a 



2 Qn 1 

-iHs/n\n - 1,R) <g> |1) 



n,fl)®(|0) + |l» 



+iv Vn+1 |T+ ® s (n-xn),-R><8 |0> 

® [ <n, J2| <8> ((0| + (1|) 
- (g) (1| 



-iisV^+l (T+® s (n-xn),R\®(0\ (104) 



All we need to do now to evaluate the final state of the 
qubit is to trace over the field degrees of freedom: 



Pq.out — Tr^ (p a 



(105) 



It is straightforward to calculate this partial trace of 
Eq. I|l 04(1 over the Fock space; the result is 

r~i2 p — 27rnSl/a 

Pq , out = _^ [(|0) + |1))((0| + (1|) 

n 

+ |H 2 n|l)(l| + ^ 2 (n+l)|0)(0|] , (106) 
where [cf. Eq. fTUU|) ] 



Introducing the sums 



S = (l-e- 27r °/ a )^ 



.,— 27r7iO/ a 



Qn 



S a = (1-e 



3 — 27rnf2/ a 



(107) 



(108a) 



(108b) 



5 e = (1 - e - 2 ^/«)^ J2 [n+ ] n (108c) 

Qn 

which satisfy S + S a /2 + S e /2 = 1, Eq. ifTTOI) can be 
written in the more compact form 



1 I So + S e Sq 
Pq ' out ~ 2 I So So + S a 



(109) 



The purity Tr(p 2 out ) can now be expressed as a func- 
tion of the acceleration; a numerical plot is shown in 
Fig.H 



VI. CONCLUSIONS 

In conclusion, we have studied the effect of quan- 
tum fluctuations in the spacetime metric on the evolu- 
tion of a qubit. We considered three general-relativistic 
paradigms: flat fluctuations in a two-dimensional 
Minkowski space, mass fluctuations in the Schwarzschild 
metric for a qubit in a circular orbit, and the interac- 
tion of a qubit with coherent and squeezed gravitational 
waves. Finally, we analyzed the decoherence of an accel- 
erating qubit interacting with the thermal bath of Unruh 
particles in Minkowski spacetime. 

We found that in addition to the expected decoherence 
of the qubit, fiat fluctuations induce a phase drift propor- 
tional to the square of the magnitude of the fluctuations. 
Similarly, mass fluctuations in a black hole induce a phase 
shift in addition to the decoherence. Gravitational radi- 
ation has a peculiar interaction with the qubit, involv- 
ing periodic episodes of decoherence followed by restora- 
tion of coherence for a sinusoidal wave. These periodic 
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episodes disappear when we consider wave packets which, 
in both the coherent and squeezed states induce a small 
phase drift proportional to the strength of the field as 
well as the expected decay of off-diagonal density-matrix 
elements. The thermal Unruh radiation induces decoher- 
ence with a characteristic dependence on the magnitude 
of the acceleration. 
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